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We study the stability of solid- and supersolid (SS) phases of a three-dimensional spin- and a hardcore-Bose- 
Hubbard models on a body-centered cubic lattice. To see the quantum effects on the stability of the SS phase, 
we model the vacancies (interstitials) introduced in the solid, which are believed responsible for the appearance 
of the SS phase, by spinwave bosons and adopt the interaction between the condensed bosons as a criterion. 
A repulsive nature of the low-energy effective interaction is the necessary condition for a second-order solid- 
SS transition and, when this condition is met, normally the SS phase is expected. In calculating the effective 
interaction, we use expansions from the semiclassical- (i.e. large-.?) and the Ising limit combined with the 
ladder approximation. The impact of quantum fluctuations crucially depends on the energy of the solid phase 
and that of the superfluid phase at half filling. As an application to 4 He, we study the parameter region in the 
vicinity of the fitting parameter set given by Liu and Fisher. For this parameters set, quantum fluctuations at the 
second order in S^ 1 destabilize the solid phase, which is supposed to be stable within the mean field theory. 
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I. INTRODUCTION 

The supersolid (SS) state, which has both diagonal- and off- 
diagonal long-range order, has been investigated over the past 
five decadesi^iS. Recently, Kim and Chan suggested 4 that the 
observed non-classical rotational inertia (NCRI)^ in solid 4 He, 
might be attributed to coexisting superfluidity. This experi- 
ment sparked a renewed interest and the origin of the NCRI is 
still under debate^. 

The quantum lattice gas model (QGM) 7 -, or equivalently the 
hardcore-Bose-Hubbard model, is one of the simplest mod- 
els suited for studying the low-temperature physics of quan- 
tum solids. Since the QGM is in an exact correspondence to 
5=1/2 quantum spin models^, we can use powerful meth- 
ods developed in quantum spin systems in understanding the 
physics underlying the QGM. The QGM has been applied^ 
to study the possibility of the SS in 4 He, and later the com- 
prehensive discussion^ given by Liu and Fisher concluded, 
within the mean-field approximation (MFT), that the SS ex- 
ists in 4 He. However, recent studies on the SS in the 2D 
square lattice systems revealed that quantum fluctuations dra- 
matically change the behavior and may even suppress the SS 
which is supposed to exist within the MF T 10 ' 11 - 12 . For the 
optimal fitting parameter set obtained by Liu and Fisher for 
4 He (LF point; see ©), frustration seems to play an impor- 
tant role. Hence, interplay between quantum fluctuations and 
frustration may change the physics of the QGM of 4 He. 

Recently, Bose-Einstein condensation (BEC) of magnons 
has been observed experimentallyi^ and is now widely 
investigated 1 -. Effects of frustration on magnon BEC would 
be intriguing in their own right, as frustration may enhance 
quantum effects and even lead to such exotic condensed states 
as the SS which are hardly realized in real Bose systems. 
For instance, quite recently, Takigawa et al. reported-^ a 
persisting spin superlattice in SrCu2(B03)2 (SCBO)^ 7 - co- 
existing with (possibly) mobile magnons even beyond the 
1/8-plateau, which is reminiscent of the SS state predicted 
theoretically 1 - for SCBO. Because weak anisotropic interac- 
tions break the rotational symmetry around the externally- 
applied magnetic field (or, U(l) gauge symmetry in the QGM 



language), this phase may not a true SS phase. However, the 
discrete subgroup of the rotational symmetry can be sponta- 
neously broken^ and the observed phase might still hold a 
close relationship to the SS in its original sense. The physics 
of this phase and the realization of the 'magnon SS' in other 
compounds are also topics to be investigated more closely. 

For the clear understanding of NCRI in 4 He and the SS 
states in spin systems, it is useful to find a criterion which 
assesses the combined effect of quantum fluctuations and frus- 
tration on the stability of the SS phase. In this paper, with the 
help of spin wave expansion, we push ahead with the widely 
accepted intuitive picture 2,3 that BEC of vacancies or inter- 
stitials gives rise to the SS state to propose that the interac- 
tion among the condensed vacancies (interstitials) serves as 
a good criterion for the stability of the SS. To this end, we 
adopt the so-called dilute-Bose-gas technique— i. Normally, 
the dilute-Bose-gas approach is used only in the vicinity of 
the saturation field to obtain unbiased (asymptotically) exact 
result a 21 ' 22 ' 23,24 ' 25 , since the lack of an exact reference state 
(i.e. vacuum) on which boson excitations are defined ham- 
pers the construction of a well-defined bosonic Hamiltonian. 
To overcome this difficulty, we introduce the spin magnitude 
S and the Ising-like anisotropy as large control parameters, 
which guarantee the validity of the reference state even far 
below the saturation field, and develop a systematic expansion 
with respect to these parameters. 

The organization of the present paper is as follows. In 
Sec. lU we introduce a three-dimensional model Hamiltonian 
on a body-centered cubic (bcc) lattice (see Fig. [TJ and briefly 
review the correspondence between the spin model and the 
QGM. At the same time, we classify the ground-state phases 
within the MFT. Then, we derive a spin-wave Hamiltonian by 
using the Dyson-Maleev transformation in the solid phase. 

In Sec. [Ill] we outline the dilute-Bose gas approach used 
in investigating the SS phase around the solid phase. If the 
effective interaction among the condensed bosons is attrac- 
tive, the SS phase for low condensate density is normally 
phase-separated. Although we do not exclude the possibility 
that the SS emerges through the first-order transition from the 
solid phase, this seems unlikely from various results obtained 
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by quantum Monte-Carlo simulation a 10 ' 1 1 1 12 - 26 i 27 i 28 . Hence, in 
this paper, the SS is said to be 'unstable' ('stable') if the inter- 
action between condensed bosons is attractive (repulsive). To 
evaluate the interaction concretely, we need approximations. 
In Sec. IIV1 we study the properties of the solid and the SS 
phase by the large-5 expansion up to the second order in S 1-1 . 
At the first order, the MFT results are reproduced. We shall 
find three types of SSs, which have properties similar to those 
appearing in the 2D-square lattic e) 27 ' 29 ' 30,31 . The formulation 
of the second-order perturbation is detailed. 

In Sec. [V] we study the properties of the solid- and the 
SS phase by the Ising expansion up to the second order. Al- 
though quantum fluctuations seem to suppress the interactions 
at the first order, the boundary determining the stability of the 
SS does not shift and the magnetization process is affected 
only quantitatively by quantum fluctuations. In other words, 
the stability itself is known from the MFT if the large Ising 
anisotropy exists. To see the effect of quantum fluctuations 
on the stability of SS, we have to proceed to the second order 
calculation. 

Our main results are summarized in Sec. I VII where we 
study the stability of the solid and the SS phases focusing on 
the LF point. Readers who only want to know the main re- 
sults may skip Sec. lIVI andlVland go directly to this section. 
For the parameter set corresponding to the LF point, quan- 
tum fluctuations destabilize the solid state expected from the 
MFT (Fig. ITTb at least within the conventional second-order 
spin-wave expansion. Concerning the stability of the SS, both 
of the two second-order calculations conclude that quantum 
fluctuations only slightly change the MFT boundary of the 
SS phase, provided that the energy of the solid phase is suffi- 
ciently smaller than that of the superfluid phase at half filling 
(Fig. [331). In the vicinity of the LF point, where the above con- 
dition is not satisfied, it is suggested that the SS phase is frag- 
ile against quantum corrections or even completely smeared 
out, although the validity of both approaches is not obvious in 
this region. 

For concreteness, we restrict our discussion in this paper to 
a quantum spin model on a bcc lattice. However, our approach 
can be easily generalized to quantum spin models on other 3D 
lattices. 




FIG. 1: (color online) Three-dimensional body-centered cubic (bcc) 
lattice and interactions considered in the text. Filled circles denote 
spins connected by anisotropic (XXZ-like) exchange interactions. 
We divide the lattice into two sublattices, which are distinguished 
by the size of spheres. Each sublattice forms a simple cubic lattice. 



H. SPIN HAMILTONIAN AND QUANTUM LATTICE GAS 
MODEL 

A. Model Hamiltonian 

Let us consider the following frustrated spin Hamiltonian 
on the bcc lattice with the nearest neighbor Ising antiferro- 
magnetic (AF) interactions ( J{ > 0): 

H = E { J i ^ s; + Jttsrsf + sfsf)} 

n.n. 

+ E { J i s ?> s l> + Jit(S$Sf + Sf,S*)} (1) 

n.n.n. 

+ shJ2s?, 

i 

where the summations n.n. and n.n.n. are taken for the 
nearest-neighbor and the second-nearest-neighbor pairs, re- 
spectively. This Hamiltonian in the case that J? = for 
i = {1,2} (Heisenberg case) has been investigated from the 
various approaches.— In the case of S = 1/2, this Hamilto- 
nian is equivalent to the following hard-core bosonic Hubbard 
model* 3 - 3 - 

H = E {i~(Pi« +PiP h + J i^j} 

n.n. ^ ' 

+ E {^(pl'Py + JZnvnA (2) 

n.n.n. ^ ^ 
i 

where nj = p\pi- This model can be used to study the low- 
energy physics of 4 He if we approximate the Bose gas by the 
QGM 7 . Specifically, the 'longitudinal' couplings Jf 2 and 
'transverse' ones 2 mimic the interaction potentials and 
the kinetic energy of Helium, respectively, and the external 
magnetic field h (or \Xh) controls the pressure. In the QGM, 
J^ 2 < and Jj 1 / 'J 2 is fixed at 1/2 because of the lattice 
structure. Liu and Fisher suggested several sets of fitting pa- 
rameters appropriate for 4 He and concluded that the stability 
of the SS phase is ensured within the MFT 9 . However, the 
existence of quantum fluctuations and frustration effects may 
destroy the classical ground state. To see the validity of the 
MFT, in Sec. [VI] we shall study these effects on the ground 
state in the vicinity of the following parameter set (LF point; 
the case (a) in Ref.|H): 

J{ = 2.60, J 2 Z = 1.59, = -1, j£ = -0.5 . (3) 

Since the Ising-like Neel antiferromagnetic (NAF) phase is 
identified with a solid phase of 4 He, we restrict ourselves 
only to the case that NAF order along the ^-direction appears 
around h = and will not consider the Ising-like collinear 
antiferromagnetic (CAF) phase which realizes, in the classi- 
cal case, when 2Jf < 3J|. Let us briefly discuss possible 
classical phases at h = 0. The classical phases fall into three 
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fundamental classes (NAF, CAF, FM) as is shown in Fig. [2] 
These phases are further classified by whether the spins align 
along the z-axis or in the xy-plane. In the former case, the 
ground state may be gapped. In the latter case, the sponta- 
neous symmetry breaking of the rotational symmetry around 
the z-axis (U(l)) induces the gapless Goldstone mode and the 
phase is viewed as a superfluid (SF). The J7(l)-broken phases 
accompanied by translation-symmetry breaking in the diag- 
onal channel (i.e. (p\p\) or (S*)) as well are thought of as 
spin-analogues of SSsi22i2ii 

When the spins align in z-axis, the energy of each Ising-like 
phase is given by 
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FIG. 3: (color online) Classical phase diagram for h = 0, Jt = 
— 1, J ^ = —1/2. NAF and CAF are implied as the Ising-like gaped 
ones. The dot labeled as LF represents the LF point ([3}- 



The magnetization curve for the LF point ([3]) is shown in 

Fig. a 



where N is the number of sites of each sublattice. When the 
spins align in xy-plane, the phases are viewed as SFs and the 
energy of each phase is given by 
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FIG. 2: (color online) Spin configurations for the three phases ('FM', 
'NAF' and 'CAF') in the text, (i) 'FM' (ferromagnetic phase) repre- 
sents a phase where all spins are polarized along the field direction, 
(ii) In 'NAF', the spins on each sublattice align ferromagnetically 
while those on different sublattices are anti-parallel, (iii) 'CAF' is 
made up of two antiferromagnetically-ordered sublattices, which, as 
a whole, align in a collinear manner. 



The ground-state phase diagram of the QGM for h 
Ji 1^2 = 1/2 is shown in Fig. [3] 



Oand 



FIG. 4: (color online) The magnetization curve obtained for the 
LF point (O within the MFT. Magnetization M is given by 
(1/2N) ^2i N (Si). The existence of the phases solid, SS and SF 
is confirmed. 



To see the physics of the NAF phase more clearly, we divide 
the whole lattice into two sublattices A and B each of which 
forms a cubic lattice. Since we can change the sign of the 
(n.n.) transverse coupling Jj 1 — > — Jj 1 at will by making spin 
rotation (by it) around the z-axis St — ► —SHi = x,y,l G A) 
only for the A-sublattice, we may restrict our consideration to 
the case < 0. 

The correspondence between the phases in the quantum 
lattice-gas formulation and the ones in the (quantum) spin- 
model formulation is shown in TABLE U In TABLE HJ the 
long-distance limit |i — j| — > oo is implied. In the supersolid 
(SS) phase and its spin counterpart, translation symmetry is 
spontaneously broken (i.e. (S? ) and (ni) modulate in space 
with non-trivial periods) as well as the axial U(l). 

In this paper, we reserve the terminology 'NAF' for the 
Ising-like NAF phase and the corresponding phase in the 
Bose-Hubbard model (O will be called (checkerboard) 'solid' 
or 'half-filled solid'. 
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TABLE I: Correspondence between the quantum lattice-gas model 
(QGM) and the spin model. 'TS' denotes translational symmetry. 



where 



QGM(Bose-Hubbard model) 



Spin model 



vacuum 
checkerboard solid 
striped solid 
SF«PiPj) + withTS) 
SS ((pipj) ^ with broken TS) 



(polarized) FM 
(Ising-like) NAF 
(Ising-like) CAF 
(S^Sf) ± withTS 
(Si+Sj 1- ) / with broken TS 



B. Spin wave Hamiltonian 

In order to rewrite the spin operators in terms of bosons, it is 
convenient to define magnons over a presumed reference state. 
In the case of NAF, all spins on the A-sublattice point upward 
and those on the B-sublattice downward (see Figf2]i. There- 
fore, it would be reasonable to introduce the following anti- 
ferromagnetic Dyson-Maleev (ADM) transformatio n 34,35 ' 36 , 



5+ = V2Sa h Sr = V2Sa\ 1 - 



~2S 



s„ 



Si = S - a\a h for Is A. 

: = V2SbL s m = V2s(i 



bLbr, 



25 



5™ S ' 



j6 m , for m e B 



(6a) 



(6b) 



If we introduce the Fourier transformation as 



a\ = 



k in 



(7) 



(N is the number of sites of each sublattice), then the Hamil- 
tonian is given by 



H =H +H X + const , (8a) 
F " = £ S {( £ o(k) - h)a{a u + (e (k) + h)^ 

k 

+ to(k)(o k 6_ k + 4bL u )} , (8b) 
H ^=^ E {- 2J i^i( c l)4 1+q 4 2 - q «k 1 6k 2 

ki ,k2,q 

+(J|C 2 (q) - J 2 ± C 2 (k 2 ))aj ii+q aj C2 _ q a kl a k2 

+ (J|C 2 (q) - #C a (k a - q))fo kl+q 6 k2 - q ^6k 2 

-J^d(k 2 ) (4 1+q a k2 _ q a kl 6i k2 +^ 1+q &k a +q6k 1 a_i (i ,)} 

(8c) 



e (k) = 8Jr - 6J| + 2J 2 - L C 2 (k) , 
t (k) = 2J 1 i C 1 (k), 

— cos -f cos — , 
2 2 2 ' 

C 2 (k) = cos k x + cos fc a + cos k z . 



C\ (k) = 4 cos — - cos -i 



(9a) 
(9b) 

(9c) 
(9d) 



Although this Hamiltonian is not hermitian and contains un- 
physical states^, we believe that the Hamiltonian given by 
eqs.(l8all-(l8c]i correctly captures the low energy physics at 
and around the half-filled solid. Actually, in the case of 
magnon BEC just below the saturation field, though gener- 
ally not proven, it is known for some specific models that the 
ferromagnetic Dyson-Maleev transformation, the Holstein- 
Primakoff transformation and the hard-core boson expansion 
for the 5=1/2 case give the same ground state in a dilute 
bose gas approac h 21 ' 22,23 - 24 . 



HI. GENERAL FORMALISM 

In this section, we outline the dilute-Bose gas approach by 
which we shall investigate the SS phase around the half-filled 
solid appearing in the system described by the Hamiltonian 
eq.®. 



A. Bogoliubov transformation 

In the following analysis, we frequently deal with Hamilto- 
nians of the following form: 



ad = 5{(e(k) - h)ata k + (e(k) + h)b\j}^ 
+ t(k)(a k 6- k + a k &l k )} . 



(10) 



This is the most general quadratic Hamiltonian allowed by 
hermiticity and sublattice symmetry. When we consider the 
quadratic part d8bb of the Hamiltonian H, the functions e(k) 
and t(k) should be taken as: 



e(k)=e (k), t(k)=fo(k) 



(11) 



However, since the interaction Hi shifts the grounds state, the 
renormalized quadratic Hamiltonian which leads to the exact 
Green's function including the self-energy do not in general 
coincide with ( l8bl . Generically the functions e(k) and t(k) 
are given by: 

e(k) = e (k) + e'(k), i(k) = t (k) + t'(k) , (12) 

where e'(k) and t'(k) are of the order of 5 _1 since the in- 
teraction Hi is of the order of 5°. In this paper, we approx- 
imately calculate the functions e'(k) and t'(k) in powers of 
5 _1 (Sec. HVI i or of the Ising coupling constant 1/ Jf (Sec.lVl). 
Now let us assume that we have found an appropriate iJ qU ad- 
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Then, in order to eliminate the off-diagonal terms ab + aJb\ 
we may introduce the following Bogoliubov transformation: 

a k = cosh^kQ^k — sinh^k/^k > (13a) 
6k = — sinh 6*k«k + cosh0k/3 k , (13b) 

which transforms iJ qua d to: 



Hqu&d — S 



{(e a (k) - h)a{a k + (e (k) + h)^ 
+ /(k,0 k )K/3-k + « k / 3 ik)} • 



(14) 



In the above, we have introduced two functions 

e a (k) = e^(k) = e(k) cosh26» k - i(k)sinh20 k , (15a) 
/(k, k ) = -e(k) sinh20 k + t(k) cosh26» k (15b) 

If we choose 6* k in such a way that /(k, 6*k) = 0, i.e. 

t(k) 



tanh 26>k 



e(k) 



(16) 



ffquad is diagonalized and reads 
H quad = S(e Q (k) - /i)4« k + S{e a (k) + h)plf3 k , (17) 

It is important to note that the magnetic field h has different 
signs for a and j3. Assuming the (unique) minimum of the 
spinwave excitation e Q (k) takes place at k = Q, we may 
introduce the renormalized chemical potential by 



p a = h- e Q (Q) 



(18) 



Now suppose we increase the external magnetic field h (or 
p a ). Then, the gap of the a ((3) boson decreases (increase) 
and the a bosons into an a-SF phase discussed below. 



B. Supersolid from magnon-BEC 

In the previous subsection, we have seen that, as the exter- 
nal magnetic field is increased, the a magnon condenses at 
fi a = while the other remains gapped. Now, we show that 
this BEC of the Bogoliubov-transformed magnons generally 
leads to an SS phase. When a BEC occurs for p a > 0, «q 
takes a finite expectation value («q) ^ and, correspond- 
ingly, the original bosons a, b have the following expectation 
values^ 

(a Q ) = cosh6>Q(aQ), (6q) = -sinh6»Q(a^) . (19) 



In a dilute-gas limit, when translated into the spin language, 
this implies the following spin configuration^ 



^coshflq cos(Q • 1 + ip){l + ^%^), 

sin(Q. 1 + + 
{S - AS) - pcosh 2 6 Q , for 1 e A , (20a) 
- sj2Spsinh6 Q cos(Q • m + tp)(l + ^^-), 

iJ 

=T\^smhe Q sin(Q-m + ^)(l + ^^), 

smh 2 6» Q , for m e B , (20b) 



(Si 

(Sf) = ±v/2S>cosh6»Q 

3 cosh 2 9q 



{Si 
(A 



(S^) = -(S-AS) + p 

where the real-space wavefunction is given by (a r ) = 
^/pexp{±i(Q ■ r + ip)} and AS = l/A^J] g sinh 2 9 q . The 
function f(AS) = AS/2 + 0(1/ S) is obtained from the 
Holstein-Primakoff transformed operator S^ and is indepen- 
dent of pq in the dilute-gas limit. One can easily see that this 
state may be thought of as an SS of magnons; an off-diagonal 
long-range (incommensurate) xy-order (which translates into 
an SF long-range order) and a diagonal (commensurate) 2- 
sublattice z-order coexist with each other. In general, a mod- 
ulation in the transverse component S x ' v with the wave vector 
Q is incommensurate with the pattern of the z-order. 

If we denote the effective two-body interaction among the 
condensed bosons evaluated at p a = 0~ by T, the leading 
term of the system energy is in general written, as a function 
of the condensate density p, as 



E e f( 



const 



Sfi a p 



(21) 



(22) 



Then, provided r > 0, p is given by minimizing E: 
a = -p for ^ Q > . 

However, the condition T > is not sufficient condition for 
the stability of the SS phase since there may be higher order 
terms with negative coefficients in E e ff, which may select a 
very large value of p and eventually destabilize the SS phase. 
If T < 0, on the other hand, one may expect a phase separa- 
tion accompanied by magnetization jump near p a = 0. For 
both cases, there exists an additional possibility of more exotic 
phases where single-particle BECs are no longer relevant^ 

The low-energy excitation spectrum of the SS phase is eas- 
ily obtained as in the ordinary superfluid Bose gas^. Defining 
k = q—Q, we may expand e Q (q) = e^+kikj / (2rrii 3 ■) + ■ ■ ■, 
where the summation over repeated indices is implied. We can 
diagonalize my to obtain a standard dispersion kikj/2m,ij = 
k'/ 2 / (2m 1 /) = e g (k'). Using this notation, the excitation spec- 
trum of the SS phase is given by 



ftss(k) = ^(k) 2 + 2^ R e g (k) a ^/25/i R e g (k). (23) 

For finite temperature, the Bose condensed bosons are sup- 
pressed, and the critical temperature is given by 

k B T c = 2.087(m x m y m z )-i(-^)i . (24) 
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For T > T c , the long-range order disappears and (5 ± ) = 0. 

Above discussions assume the dilute-gas limit, where the 
scattering length is much smaller than the average inter- 
Specifically, our approximation is 



atomic distance 
valid when 



-1/3 



T(m x m y m z p) 1/3 < 1 



(25) 



To summarize, the knowledge about the wave number Q at 
which the magnon BEC occurs, the effective mass m, and the 
effective (2-body) interaction T for the condensed bosons en- 
ables us to derive the stability, the spin configuration which is 
not commensurate with the assumed sublattice structure, the 
quasi-particle excitation spectrum and the critical temperature 
of the SS phase. Therefore, the analysis boils down to the 
calculation of Q and T. A remark is in order here about the 
definition of the bosonic vacuum. In eq.©, it is implicitly as- 
sumed that the NAF phase gives a well-defined vacuum (i.e. 
the ground state when the condensate is absent) for the two 
bosons. In general, the NAF state shown in Figffjsuffers from 
quantum fluctuations and the above assumption is justified ei- 
ther for the semiclassical (i.e. large-S*) case or the Ising-like 
(i.e. large- Jf/J^) limiti! 

In the following sections, we carry out the calculation by 
combining the ladder approximation with the large-S* and the 
Ising expansions. Concretely, in Sec. IIVI we will obtain 
eq.([T2b and the interaction T by the large-S* expansion up to 
the second order in S^ 1 . At the first order, our approach will 
reproduce the results of the MFT; there are three types of SSs. 
At the second-order perturbation, quantum fluctuations may 
change the properties of the solid and the SSs qualitatively. 
However, we will see that the large-S" expansion is not reli- 
able to calculate T when the Ising-like anisotropy Jf is large. 
To overcome this difficulty, we will study eq.([T2b and T by the 
Ising expansion up to the second order in Sec.[V] At the first 
order, quantum fluctuations suppress the interactions, but the 
stability of the SS itself is known by the MFT. Although we 
will find the stable bound-magnon state, this condensed phase 
may be phase-separated for large Jf . In the second order, we 
will see the effect of quantum fluctuations on the stability of 
SS clearly. 



the quadratic part of the Hamiltonian is diagonalized up to 

O(S) (see eq.fO). We note that 9 k is well-defined when 

| tanh26'^ 1) | < 1. Concretely, the half-filled solid is stable at 
h = when 



4 Jf - 3J| + 3 Jt >4|Jf 



(27) 



If this inequality is not satisfied, the spins align in the xy- 
plane (SF) (see eqs.© and ©). Meanwhile, even when the 
classical ground state is CAF (2 Jf < 3 J|), this inequality 
may be satisfied and then the metastable NAF phase against 
the one magnon fluctuation may be obtained. In this paper, 
we will not discuss the CAF case any more. 

Let us discuss the minimum of the dispersion e Q (k) = 
e^(k) (see eq. dl5al )) to determine the structure of the SS. 
From 1261 . the dispersion relation reads, 



Ji 



V(k) = ^o(k) 2 -io(k)^ 



(8Jf - 6Jf + 2J 2 - L C 2 (k)) 2 - 4J 1 L2 Ci(k) 2 . 

(28) 



where the superscript (i) of e^k) denotes that the function 

is evaluated at Ok = 6 k . In the following we shall use this 
notation to the other arbitral functions of 6>k- The minimum is 
obtained by setting Q = = (0,0,0) or Q 2 = (tt,tt,tt). 
Although we can not exclude other possibilities generally, this 
is always the case for the parameter sets used in this paper. It 
is convenient to introduce A as 



A = 4i ) (Q 2 ) 2 -4 1 i ) (Qi) 2 = 16(-J 2 x (12J?-9Jf)- 



Ut 2 ) , 
(29) 

Then, one chooses k = Qi when A > or k = Q 2 when 
A < 0. We have checked that the SS with Qi (SSI) is always 
favored for < 0. And when the Ising-like anisotropy is 
large, i.e., 12 — 9 Jf ^> 4Jf~, very small positive selects 
the SS with Q 2 (SS2). For each case, we plot the dispersion 
relation esi(k) along the (1, 1, 1) -direction in FigfS] 



IV. PERTURBATION THEORY IN S 

In this section, we study the physics of the SS phase by 
the perturbation theory in the parameter S . The first-order 
calculation gives the same ground-state phases as the MFT. At 
the second order, on the other hand, quantum fluctuations play 
an important role and may destroy the classically stable solid 
(NAF) or the SS phase. 




A. First-order perturbation 

If we assume 6*k by 
*o(k) 



tanh26»^ 1) = 



e (k) 4Jf-3Jf + J 2 1 C 2 (k) 



(26) 



FIG. 5: (color online) The dispersion relation of the excitation energy 
esi(k) for k = (k, k, k), h = 0, Jf = 3, Jf = 1, Jt = -1. G) 
is obtained at =0.1 and the minimum is at k = Qi. (ii) is at 
J-2 = 0.148 and the minima are at both k = Qi and Q2. (iii) is at 
J 2 =0.2 and the minimum is at k = Q2. 
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From eq.dT8b. the chemical potentials of both phases are 
given by 



MsiSSi = h - esi(Qj) , fori = {1,2} 



(30) 



where the subscript 'Sn means that the interactions are ex- 
panded up to n-th order in S^ 1 and 'SSi' represents the types 
of the SS. The effective masses are isotropic and are respec- 
tively given by 



y/(4Jf - 3Jf + 2,Jj) 2 - 16Jp 
TOS1SS1 = —77 TTTT-T, „ T , , 7m~\ . . T i9i > ( 31a ) 



2S{-J 2 - L (4Jf - 3J| + 3^) + 4^2} 
1 



™S1SS2 



2SJ 2 L 



When A = 0, the two minima are degenerate and we have 
to take into account two independent condensates and phases 
which are not characterized by (fTTI i may appear. A brief dis- 
cussion on this case is given in appendix lAl The SS phase of 
4-sublattice structure (SS3) actually exists for certain param- 
eter sets. There exist three types of SS around the half-filled 
solid. 

Next, we consider the stability of the SS phase. The 2- 
body interaction between a bosons is given by the first-order 
diagram since the bare Green's function of a {(J) bosons is 

i/{io — Sea (k) ± h) — 0{S^ 1 ) for u> ~ — fi a and the vertex 
function is 0(S Q ). The alternative view is that, if we rescale 
the Hamiltonian by the vertex function is 0(S~ 1 ) and 
the diagram is suppressed by S^ 1 for each vertex. Therefore, 
we need only the vertex function between a bosons. By re- 
placing ak — * cosh^l ak and — ► — sinh^^a^ in H, the 
interaction term of a bosons appears as the following form: 



^L^K( q ; kll k 2 )a ki 



+q a k 2 -q a kl"k 2 



(32) 



where the factor 2 in front of N is considered for the symmetry 
factor. For the case A > and Q = Qi(= 0), Y is given by, 

Tsissi = V a (0; Qi, Ch) = 6(J| - j£) , (33) 

Thus, SS phase of Qi is stable for J| — J 2 > 0. 
For the case A < and Q = Q2, Y is given by, 

r S iss2 = V a (0; Q 2 , Q a ) = 6(J| + J 2 X ) , (34) 

In this phase, the Spin on B-sublattice does not have the trans- 
verse magnetization even for /^siss2 > since sinh = 0. 

To see the validity of the above picture, we compare the 
above result with that of the MFT. Although there exists the 



extensive MFT calculation of this model for» 



8.y 



< 0, to 



the best of our knowledge, there is not the appropriate mean 
field calculation of the models for J 2 > 0. Hence, we redo 
the MFT for S = 1/2. Now, the ground state energy is ob- 
tained by replacing the operators in H with their expectation 
values of Pauli matrices on each site, e.g., Yluj) ^i^j ~ ¥ 
j) S 2 { a j) ■ We compare energies of the three types 



of spin configurations, 

El 



NS 2 



- = 8J?(a z )(a z )' + 3Jf((a 2 ) 2 + (a z )' 2 ) 

8| J^\tt' + ZJi(T 2 + r' 2 ) + h((a z ) + {a*)'), 

(35a) 



rrrnean 

^ = 8J z (a z )(a z Y + 3Ji({a z ) 2 + (^)' 2 ) 

^2 1 _/2 



3Jt(r A +r") + h((a z ) + (a z )') 



E 



l/4fflled 

NS 2 



-h. 



(35b) 
(35c) 



(31b) where t = ^f{(T x ) 2 + (av) 2 = y/l - (a z ) 2 . E™ ean is ob- 
tained from the 2-sublattice structure, i?™ 63 " is from the 2- 
sublattice structure of (a z ) and r with AF-(tt,tt,tt) (a x ' v ) 
ordering on each sublattice, and -Ei/4fin e d is from the quarter- 
filled solid. In this paper, we ignore the possibility that another 
types of SS phases appear around the quarter-filled solid as in 
the model on the square lattic e 12 ! 31 . By minimizing each en- 
ergy numerically, we obtain magnetization curves for various 
parameters. We confirmed that the Bose-gas approach gives 
the same results as the MFT one. The specific examples are 
shown in FIG|6] 



+ 

CO 




-h 



FIG. 6: (color online) Magnetization curves obtained from the MFT 
for Ji = 3, J\ = — 1. We assume 2-sublattice structure or the 
quarter-filled solid. M is given by (1/27V) J2t (Si) ■ (0 

curve is 

obtained for Jf = 0.5, j£ = 0.3. (ii) is for Jf = -0.1, J 2 X = 0.4. 
(iii)isfor Jf = 0, J2 = -0.4. (iv) is for Jf = -0.5, j£ = -0.4. 
All curves has a half-filled solid around h = and SF phase below 
the saturation field, (i) curve has SS2 just ahead the half-filled solid, 
and the quarter-filled solid, (ii) has SS2 also, (iii) has SSI which is 
connected to SF phase continuously, (iv) does not have SS phase. 



B. Second-order perturbation 

To see the effect of quantum fluctuation more clearly, we 
consider the second-order perturbation theory in the param- 
eter S~ l . In this order, the ground-state phase may become 
different from the mean-field one. 

To begin with, let us consider the state of the half-filled 
solid by diagonalizing the quadratic term in the Hamiltonian. 
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When the interaction terms ( I8cb are written in terms of the 
Bogoliubov-transformed bosons and put into the normal-order 
form, additional quadratic terms appear. As a result, the 
quadratic part of Hamiltonian reads, 

-ff q uadS2 ={(Se (k) - Ti(k))cosh2# k 

- (5*o (k) - T 2 (k)) sinh20 k - Shjala k 
+ {(5'eo(k) -T 1 (k))cosh26» k 

- (Sf (k) - T 2 (k)) sinh20 k + Sh^plfa 

+{-(Se o (k)-T 1 (k))sinh20 k 

+(5to(k) - T 2 (k)) cosh 20 k ) }(a k /3 u + 4/3+) , 

(36) 

where XfcS are given by eq. dBll ). Even in the normal-ordered 
2-body interaction terms, there exists the terms which shift 
the vacuum with respect to a and (3 (e.g. ofi /3* aJ /3' |0) ^ 0), 
which leads to the self energy. However, this contributes the 
Green's function in the third-order of S 1 " 1 and we neglect 
the self energy in our approximation. We note that, even if 
we use Holstein-Primakoff transformation, the same quadratic 
Hamiltonian is obtained up to the second order in 3 _1 . The 
difference between the two boson representations (i.e. Dyson- 
Maleev and Holstein-Primakoff) appears in the 2-body inter- 
action term. 

Now, 6> k is given by solving 



(-Se (k)+Ti(k))smh26> k 
+ (St (k) - T 2 (k))cosh26» k 







(37) 



To evaluate I\, 2 , we need the explicit form of the function 6* k . 
Since Ti j2 is suppressed by a factor 1 / S in the diagonalization 

procedure, we use 6^ which is obtained in the first-order cal- 
culation for the integrands in eq. dBU . Therefore, 6 k 
is corrected up to second order is given by 



(2) 



which 



e (k)-T 1 (1) (k)/5 



(38) 



(2) 

If |tanh26>n > 1, the spinwave expansion concludes that 
the half-filled solid is unstable and that other phases may take 
over. In fact, this happens for certain choices of the parame- 
ters. The detailed result will be discussed in seclVIl Then the 
quadratic Hamiltonian and the dispersion relation es 2 (k) are 
given respectively by 



S (e S2 (k) - h) alau + S (e S2 (k) + h) /3 k /3 k , (39) 



If we introduce the appropriate constants a\,...,a^, the 
above phonon dispersion es2(k) may be written generally as: 



e S2 (k) - V(ai+02C 2 (k)) 2 -(a 3 Ci(k)) 2 



(41) 



and qualitatively the same dependence on k as in the first- 
order case is obtained. In our calculations, the minimum is 
always locked at Qi = (0,0,0) or Q 2 = (tt, 7T, if), which 
respectively corresponds to SSI or SS2. The criterion, which 
determines the structure and the effective mass for each phase, 
is easily obtained in the same manner as in the first-order case 
(see eq.(|29l and (|3T1)). However, the explicit forms are some- 
what lengthy and we do not show them in this paper. In the 
following, we shall concentrate on the physics of SSI and SS2 
and shall not discuss SS3 further. The chemical potential [i a , 
which controls the onset of BEC, are also different from the 
first-order one ( f30b and is given by 



MS2SSi = h- e S2 (Qi) , for i = {1, 2} 



(42) 



Next, we briefly recapitulate the method by which we 
calculate the effective interaction T among the condensed 
bosons. We simply evaluate the diagrams up to the second 
order in We have one diagram at the first order and six 
at the second order. The second-order diagrams are shown in 
Fig- El To evaluate the second order diagram, we use the bare 
Green's function at ^s 2 ssi = in the dilute bose gas approxi- 
mation: 



(T(o*4)H) 



u - S (e S2 (k) - e S2 (Qi)) + *0+ 



(43) 



In the presence of a finite condensate |(a)| 2 = p oc [i, the 
Green's function, which is obtained for a new operator a' = 
a — (a), gets modified continuously from the one at the onset 
of BEC^. Specifically, (aW+) at M>0 = ^ +O(p) , 

and (a'a') a t ^>o = O(fi). In short, the modified quadratic 
Hamiltonian and the effective interaction T calculated above 
tell us the stability and the low-energy physics of solid and SS 
phase. The detailed results are shown in Sec. I VII 




o 



e S2 (k) 



\ 



co(k) 



S 



T 2 (1) (k) 



- to(k) - 



(40) 



FIG. 7: The second-order (one-loop) diagrams in S" 1 . Straight lines 
(wavy lines) denote ot (/3) bosons. Broken lines denote the momen- 
tum transfer at the interaction. 
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Finally, we consider the validity of the expansion of the 
exponential in powers of the interaction terms in the path inte- 
gral when we calculate T. If we were able to take into account 
an infinite number of terms, the expansion would be correct. 
However, now we sum up only a finite number of terms. Thus, 
we need a criterion, even though naive, for determining the 
reliability of the expansion. A natural candidate may be the 
magnitude of the expanded interaction terms. To see this ex- 
plicitly, we consider the following simple boson model on the 
simple cubic lattice: 



H 8 



^— ' 2m 

k 



k2 ,t J 1 V- 



2A< +q < 



ki,k2,q 



(44) 

In this model, the low-energy effective interaction T s be- 
tween the condensed bosons is exactly obtained as T s = 
2A/(1 + (2/n)m\) = 2A£ n (-(2/7r)mA) n . The dimen- 
sionless constant mA captures the magnitude of the expanded 
interaction terms. Thus, for general lattice boson models, we 
may expect that (mass) x (coupling constant) gives a simple 
criterion for the validity of the expansion. 

Let us apply the above criterion to our case. For the bo- 
son masses, we use (|3TT l, which are correct up to the first or- 
der in S 1-1 , for simplicity. An appropriate choice of the cou- 
pling constants may be Jf and for i = {1, 2}. For SS2, 
the criterion reads Jfm2 C i = Jf /(25J^). Hence, however 
large the spin S may be, the series expansion of V eventually 
diverges for relatively large Ising anisotropy. Similarly for 
SSI, the perturbation expansion is not converging for large 
Ising anisotropy since mi c \ ~ — 1/(25 J^~)- We have one 
more problem in the evaluation of T; when the energy dis- 
persion at the solid is nearly gapless (i.e., tanh2#k = o ~ 1 
in eqs.([T31l and (fT6]l), cosh 26q and sinh 29q have large val- 
ues(for tanh26> — > 1, 6 — > oo). We note that these problems 
are peculiar to the evaluation of T and the low-energy physics 
of the solid (NAF) is well understood by the large-S* expan- 
sion. 

From the above discussion, we may conclude that the SS 
phases obtained within the MFT, which do not change even 
after the first-order 1/S -correction is taken into account, might 
be destroyed at higher orders by quantum fluctuations. Since 
the perturbation expansion described above is ill-behaved for 
large Ising anisotropy, we have to take another approach to 
closely investigate the fate of the SS phases. In the next sec- 
tion, we shall introduce another perturbation theory with re- 
spect to large Ising anisotropy. A reliable treatment of T for 
the case with tanh 2#k=o ~ 1 remains to be an open problem. 



A. First-order perturbation 

If we diagonalize the bare quadratic Hamiltonian Hq (18bt . 



tanh 2#k 

Then, 

cosh 0k 



4Jf - 3 Jf + J£C 2 (Il) 



0(1/ Jl) , (45) 



l + 0((Jf)- 2 ), sinh0 k = O((J 1 2 )- 1 ). (46) 

If we assume that the exact #k obtained by eq.dToT) has the 
same property, the self-energy contribution to the quadratic 
Hamiltonian (e'(k), t'(k) in O) is up to 0((Jf )°) and the 
dependence on Jf of 6*k is maintained as d46b . Therefore, the 
leading-order Hamiltonian in Jf reads 

Hn = S(eo(k) - h)a{a^ + £ S(e (k) + h)^ , 

k k 

( J 2 2C, 2(q) - ^ 2 L C2(k2))a kl+q alc 2 -q Q;k i Q!k 2 • 

(47) 



ki,k 2 ,q 



where we neglect the 2-body interaction term containing j3 
bosons since the gap of j3 boson is 0( Jf ) when the gap of a 
boson closes. The meaning of the subscript 'In' is similar to 
that of 'Sn' in the previous section; it means that terms are 
kept up to 77-th order in the Ising expansion. The minimum of 
the dispersion is obtained at Qi = (0, 0, 0) for < (SSI) 
or Q2 = (tt, it, 7r) for > (SS2). The chemical potential 
and the effective mass are respectively given by 



6J 9 2 



/in = h - (8 J{ - u„ 2 



TOll 



6|#| 



25|J 2 J 



(48a) 
(48b) 



Next, let us evaluate the interaction T among the a bosons. 
Since both the Green's function of a bosons and the coupling 
constants of interaction are 0(( Jf )°), the all-order diagrams 
equally contribute to T, which is given by the ladder diagram 
(Fig|8]l. The ladder diagram T evaluated at the solid satisfies 



T(q;ki,k a ) = ^(qjki.ka) 

1 ^ T(q';k 1 ,k 2 )A"(q-q';k 1 



N ^ 
q' 



q, k 2 q) 



w(ki + q') + w(k 2 - q') - w(ki) - w(k 2 ) ' 



(49) 



where the particle corresponding to the external line is as- 
sumed to be a real one with the energy w(k) — fx. One obtains 
the parameter T for SS; (Qi) as Tss; = T(0, Qi, Qi). 



V. PERTURBATION THEORY IN LARGE ISING-LIKE 
ANISOTROPY 

In the limit J[ /* 00, the system behaves like the Ising 
model. In this section, we compute V by the perturbation the- 
ory in ( Jf Specifically, we develop an expansion in small 
coupling constants (Jf , , J% ). 







r*~i r*~ 


\ K j + 


T 


\ K i 


->J U- ->- 







FIG. 8: Ladder diagram: T represents the ladder diagram and K rep- 
resents the kernel which is not reducible to the product of 2-perticle 
Green's function. 
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Now, the kernel K and the energy u are given by 
#ii(q;k i; k 2 ) = 2JfC 2 (q) - Jj-(C 2 (ki) +C 2 (k 2 )), 



where 



wn(k) = 2S(J 2 L C 2 (k) + 3| J 2 J 



(50a) 
(50b) 



The self-consistent equation (|49i l for the ladder diagram is ex- 
actly solvabl e) 21 ' 22 and we obtain 



ss, 



6(Jf 



l#l 



0.258 



(J£+\Jl\) 
S\Jf\ 



(51) 



for i = {1, 2}. If the limit S — ► oo is taken, T reduces to the 
first-order result in the large-S expansion C( f33T > and d34b). For 
finite S, T is suppressed by quantum fluctuations. However, 
concerning the stability of the SS phases, there is no difference 
from the MFT result as far as the denominator is positive. We 
note that if 



1+0.258 



(Jj + \Jj-\) 

s\Ji-\ 



= o, 



(52) 



the effective interaction Tnss; diverges. This suggests the 
possibility of the SS accompanied by the bound-magnon 
BEC. A brief discussion on this issue will be given in 
Sec.lVCl 



B. Second-order perturbation 

In this section, we study the SS phase around the solid in 
the second-order perturbation in ( Jf ) ~ 1 with the help of large 
S expansion. 

To obtain the renormalized quadratic Hamiltonian, we per- 
form the Bogoliubov transformation and normal-order the in- 
teraction term, sorting out the terms based on d46b . Since 
near the boundary of NAF-CAF transition quantum fluctu- 
ation may play an important role, we keep only terms of 
order 0(J|/(Jf) 2 ) (we neglect the order 0((Ji/Jf) n /Jf) 
terms with n > 2). Even on the classical boundary, 6 Jf = 
(l/2)8Jf and the Jf times coordination number is suppressed 
by the large-anisotropy Jf . Therefore, the expansion may 
work. Since the renormalized 6*k satisfies eq.(l46b. the off- 
diagonal part of quadratic Hamiltonian is given by 

%2)off = 5Z{ 5 (~ e °( k ) sinh20 k + *b(k) cosh20 k ) 

k 

<jf E c i fa) sinh W ( k ) } + 4/3-k) 



i 

2 



(53) 



Now, approximately sinh #k = to(k)/2eo(k) + Ak, where 
Ak is 0(Jf _1 ). Then, the leading order of Ak is obtained 
and, as a result, 9^ is given by 



sinh0 k I2j =Aid(k) 



cosh6>[ I2) = 1 



A 2 



(54a) 
(54b) 



A, 



2J i L 5 



45(4 Jf - 3Jf ) - Jf 



(55) 



Now, the quadratic part of the Hamiltonian reads 



#(i2)o = J>i2(k) - h)ala k + (e I2 (k) + h)^ , (56) 



where 

eia(k) = 85" Jf - 65* Jf - 16 Jf A\ + 12 Jf Aj 

+ (2SJ 2 ± + 2Jfv4 I 2 )C 2 (k) 

+(2(85Jf - 65Jf - Jf)A\ - 4S'J 1 L Ai)Ci(k) 2 . 

(57) 



Then, the minimum of dispersion is obtained at Qi = (0, 0, 0) 
or Q 2 = (it, it, 7t) as in the first-order result. The chemical 
potentials and the effective masses for each phase are given 
by the same way as in Sec. lIV AR see eq.j30l> and J3TT)). 

Let us evaluate the 2-body interaction V between the con- 
densed bosons. As in the first order case in 1/ Jf , we need to 
calculate the kernel in the ladder-diagram (see FiglU. When 
the gap of the a boson closes, that of the /3 bosons is 0( Jf ). 
Then, the correlation of the (3 boson remains short-ranged for 
low energies and {T{fi^){E w 0)) = 0(1/ Jf ). The effect 
of (3 operator in the interaction term is at most 0(l/y / Jf). 
Therefore, the interaction part of the Hamiltonian which affect 
the kernel is obtained and is given in (IB2b . Now, we shall eval- 
uate the kernel. Before doing so, a remark is in order; at the 
second order of 1/ Jf an infinite number of diagrams appear 
in the kernel. Hence, with the help of large S expansion, we 
keep the term of the third-order of S (up to 2-loop diagrams) 
and neglect the term of 0{{ Jf ) _1 5 -3 ). In the selection of di- 
agrams which contribute to the kernel, we do not view Jf as a 
special contrary to the case of the quadratic Hamiltonian, for 
simplicity. As a result, in the second order of S, four diagrams 
and, in the third order of S, fourteen diagrams contribute to 
the kernel. The one-loop diagrams are given by 3 ~ 6 shown 
in Fig|7] and the part of the 2-loop diagrams are shown in 
Fig|9] When we evaluate the diagrams, we drop the term of 
0((Jf ) -2 ) after the frequency of the propagator is integrated 
out. In this calculation, we maintain the terms Jf which are 
readily obtained in the quadratic Hamiltonian, even if the con- 
tribution of these term is 0(Jf (Jf )~ 2 )- Concretely, we use 
the gap of (3 boson as 25(8 Jf — 6 Jf ) and maintain Jf of OBT ) 
and d57| i. We solve d49b by substituting the obtained kernel, 
and the interaction T between condensed bosons is obtained. 
The detailed results are shown in Sec. lVIl 
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+ 




FIG. 9: The part of 2 loop diagrams which contributes to the kernel 
in the order in 5 ,_2 (Jf ) _ . Straight lines (wavy lines) denote a (/?) 
boson. Broken lines denote the momentum transfer at the interaction. 



C. Possibility of bound-magnon BEC 

We briefly comment on the possibility of stable bound 
states. From the viewpoint of the Bethe-Salpeter equation, 
the two-particle Green's function contains the ladder diagram, 
and the divergence of T implies the existence of stable bound 
states. In fact, this method has been successfully applie d 23,24 
to search for the stable bound-magnon state in the vicinity of 
the saturated ferromagnetic phase of ID- and 3D frustrated 
magnets. Hence, the above ladder-approximation may be also 
applied to study the bound-magnon BEC around the solid. 
Now, we note that within the first order perturbation the ef- 
fective interaction Fuss (eq.dSTt) diverges when 



1 + 0.258 



(Ji 



\Ji-\ 



S\Jj 



= o 



(58) 



In fact, we found that T had a pole below the two-particle 
threshold when the left-hand side of eq.(|58]) is negative. 
Therefore, when the denominator of ( BTT l is equal to or 
negative, we may expect that, instead of the usual magnons, 
the bound magnon condenses at the center-of-mass momen- 
tum K = provided that the chemical potential is prop- 
erly tuned. We found that the energy of the bound state at 
K = (ir,ir,ir) is higher in energy than the one at K = 
and will not affect the critical value of ijl (or h) at which the 
bound-magnon BEC occurs. However, the extensive studies 
on hardcore bose Hubbard models by quantum Monte-Cairo 
simulation s 10 ' 1 1 ' 1 2 i 26 - 27 i 28 have never indicated the existence of 
an SS accompanied by a bound-magnon BEC. Hence, it would 
be useful to reconsider this problem from the energetic point 
of view, even though rough. 

Of course, some other phases may compete with the bound 
magnon phase. In particular, the SF phase, which may appear 
from the solid phase via the first-order transition, would be 
an important candidate. In the MFT, a solid-superfluid (S-SF) 
transition occurs at 



tSSF 



= 2 



16Jf - (3J| 

8J? ~ locjr 1 )! 



4|J^|-3J 2 J 



(59) 



From (148al i. on the other hand, one sees that the bound 
magnon BEC starts at h b « 8Jf - 6(J| + | = 8Jf + 
0((Jf )°) when eq.d58l is satisfied. In the case of the attrac- 
tion (Jf + \J^~\) < 0, one sees that hb > ft-s-SF and that a 
direct first-order S-SF transition occurs before the condensa- 
tion of bound magnons. Hence, although an exotic SS phase 
brought about by the bound-magnon BEC may be expected 
(note that the gap of a bound magnon closes earlier than that 
of a single magnon) in the vicinity of the solid phase, what we 
actually have is a phase separation. 

Therefore, in order to see the bound-magnon BEC around 
the half-filled solid, it may be necessary that higher-order 
terms in the perturbation in 1/ Jf shift the critical value /i by 
0((Jf )°). In this case, the approximation used in this sec- 
tion is beyond the scope of application to search the bound- 
magnon BEC. We do not go into more detailed discussion 
about the bound-magnon BEC in this paper. 



VI. PHASE DIAGRAM 

In Sec. [IV] and [V] we have described the two methods of 
calculating the minimum Q of the dispersion by which the 
spin structure of the SS phase has been determined. On top 
of it, the mass, the chemical potential, and the interaction V 
which determines the stability of the SS have been computed. 
In this section, we show the detailed results on the phase dia- 
gram paying particular attention to the parameter set (0. 



A. Stability of a half-filled solid a spinwave analysis 

The solid phase is stable when the energy gap is finite. 
Quantum fluctuations shift the energy gap and, in certain 
cases, the gap may close. In this subsection, we study the 
properties of the half-filled solid by the conventional spin- 
wave theory up to the second order in S" -1 . Even for S = 1/2, 
the approximation may work since the ground state is ordered. 

At the first order in S^ 1 , the energy gap closes even at 
h = when | tanh#£2. | = 1, and then the energy of the 
solid and the SF phase (or, a phase with magnetic long-range 
order in the .xy-plane) is degenerate within the MFT (see (0]), 
© and d26)). For | tanhflj^J > 1, the SF phase is stabilized. 
In the second order in 5 _1 , the quantum fluctuation shift 6ji 

(2) 

to 6^ and the boundary where the gap closes also changes. If 
the transition is a usual second-order one, the emergent phase 
may be SS. The first-order transition to the SF near the bound- 
ary may be also expected. However, in the case of 2D-square 
lattice, the quantum Monte-Cairo simulations indicate that at 
the Mott-SF transition point, SU(2) symmetry dramatically 
restores^ as in the classical case. Even though there exists 
the difference of the dimensionality, we may not exclude the 
possibility that on the phase boundary SU(2) symmetry re- 
stores. 

To see the properties of the resultant phases more clearly, 
we carry out the Holstein-Primakoff transformation starting 
from the SF phase (the xy-ordered NAF phase in the spin 
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language) for > and calculate the magnon dispersion 
relation in the SF up to the second order in S^ 1 at h = 0. 

Let us briefly discuss some technical aspects of the calcu- 
lation. There are two types of excitations: one is the gapless 
Goldston mode and the other is a massive (gapped) mode. To 
obtain these, we need to integrate out the functions of 9^ as 
in (IB 11 1. If we substitute the 9^ obtained in the first order as 
in Sec. lIV Bl the gapless Goldston mode remains gapless. As 
is well known, the spin-wave expansion is well-behaved if the 
ground state is classically stable. Hence, we do not extend 
the calculation to the region where the corresponding phase is 
unstable in the MFT. Concerning the gapped mode, when the 
gap closes at the first order in the solid and the SF are de- 
generate in energy within the MFT. Since the gap of this mode 
is affected by the quantum fluctuations, the phase boundary is 
shifted in the second order in S -1 . The resulting phase may 
be either the solid or the SS phase. 

As a result, a shift of the phase boundary is found in each 
phase, as is shown in Fig(lO] and FigQl] As has been dis- 
cussed above, since the spin-wave expansion is well-behaved 
in the case that the selected phase is the classical ground state, 
only on the classical phase boundary, we can compare the dis- 
persion of each phase explicitly. At the first order in both 
dispersions are gapless. At the second order in S" -1 , when 
the dispersion of one phase (solid or SF) is ill-defined, that 
of the other phase obtains the finite gap. If the system has 
the global SU(2) (rotation) symmetry, the dispersion remains 
gapless. As a result, the shifted boundary forms the almost 
straight line which intersects that of the MFT at the parameter 
set where SU(2) symmetry exists. 




0.2 0.4 0.6 0.8 



FIG. 10: (color online) Phase boundary between the half-filled solid 
(or, the Ising-like NAF phase) and the SF phase (or, a phase with 
magnetic long-range order in the sy-plane) obtained by the large- 
S expansion up to the order of S°. Te values 5* = 1/2, jf = 
1, J2 =0.5 and h = (half-filled) are used. The solid line (blue) 
for Jf > 0.5 denotes the boundary where the gap of the solid phase 
closes. The solid line (red) for Jf < 0.5 denotes the boundary where 
the gap of the massive mode of the SF phase closes. The broken line 
(black) denotes the classical boundary between NAF, CAF and SF 
phases. The dot represents the Jf = 1 and Jf = 0.5, where the 
system has SU(2) symmetry and the boundaries intersect. In the 
highlighted region (green), the dispersion has a non-zero imaginary 
part and is ill-defined. For Jf > 0.5, the emergent phase may be 
either the SF or the SS. For Jf < 0.5, the emergent phase may be 
either the solid or the SS. 




FIG. 11: (color online) Phase boundary obtained by the large-5 ex- 
pansion up to the order of S° between the half-filled solid and the 
SF phase for S = 1/2, jf = -1, Jf = -0.5 and h = (half- 
filled). The solid line (blue) denotes the boundary where the gap of 
the solid phase closes. The broken line (black) denotes the classi- 
cal boundary between NAF, CAF and SF phases. The dot labeled 
as LF (blue) represents the LF point y), which is suggested for the 
fitting parameters^ of 4 He. In the highlighted region (green), the 
spin-wave expansion is ill-defined and the emergent phase may be 
either the SF or the SS. The straight line and the broken line intersect 
at Jf = 1, Jf = —0.5, where SU(2) symmetry exists. 



Next, let us discuss the application to 4 He as the QGM. As 
shown in FigJTT] at the LF point, the solid phase is unstable 
even at h — and the resulting phase may be either the SF or 
the SS. To conjecture this phase, we plot the tanh 6>k=o of the 
massive modes as shown in FigQ~2] If | tanh6*k=o| = 1, the 
energy gap closes. With the help of the fitting line, we see that 
on the parameters (3), the SF phase may be stabilized. Hence, 
on this parameter set, the QGM does not make a sense, and 
the fitting parameters for 4 He must be reconsidered by taking 
into account the quantum fluctuation. 

l-|tanht)| 
0.5 



0.4 
0.3 
0.2 
0.1 




-LF 
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FIG. 12: (color online) 1 



I tanh 8^l | obtained of the order of S° 



for S = 1/2, Jf = 1.59, Jf = ±1, 
function of Jf If 1 - tanhfl^J = 



J 2 = — 0.5 plotted as a 
the gap closes. The solid 
line (red) is given by the massive mode on the SF. We obtained it 
for Jf = 1, which transforms to Jf = — 1 by the gauge trans- 
formation. The dashed line (blue) is obtained on the solid phase. 
The vertical line labeled as LF represents the LF point y). On the 
classical boundary, the gap of the SF phase largely opens. If we in- 
troduce the fitting line (the non-labeled broken line), the gap of the 
SF phase seems to be open on the LF point and the stability of the 
SF is implied. Moreover, the gap seems to be maintained over the 
point where the gap of the solid phase closes. Hence, the solid-SF 
first-order transition is expected at 2.7 < Jf < 3.2. 
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B. stability of supersolid 



The fitting parameters for 4 He shall shift from ([3]). Al- 
though the shift may be quantitatively large to see Fig [12] that 
is still expected to be perturbative since the quantum fluctu- 
ation is treated as a perturbation. Then, since the effective 
interaction T obtained within the first-order in S^ 1 (or the 
MFT) on (O is robust (see eq.d33l). the perturbative shift of 
the fitting parameters shall not affect the stability of the SS 
within the MFT. Therefore, we study the quantum effect to 
the stability of the SS near the LF point (0 by the T obtained 
in the second order perturbation in S^ 1 and 1/Jf. The re- 
pulsive nature of the effective interaction T(> 0) suggests the 
stability of the SS phase. 

Since J^/J^ is fixed at 1/2 in the QGM, we plot T as a 
function of Jf and Jf, as is shown in FigfOland Fig[14] On 
the phase boundary, the perturbation theory in S 1-1 gives the 
divergence to — oo because of for tanh0k=o —* 1 an d the 
used approximation is beyond control. Near ([3]), the pertur- 
bation theory in 1/ Jf also has the problem of accuracy since 
the suppression of the expansion parameter 1 / Jf may not be 
sufficient. However, both methods lead to the one identical 
conclusion. In the case of S = 1/2, both predict that T is 
considerably suppressed near (0, and, as a result, the second 
order term has the same magnitude as the relatively large first- 
order term. Hence, quantitatively, it may be understood that 
the r of the MFT (that of the first order in S^ 1 ) is not reliable 
near (01 and there exists the possibility that the quantum fluc- 
tuation breaks the stability of the SS. Therefore, even if the 
shift of the parameter set from (O is perturbative, the stability 
of the SS phase of 4 He remains to be a question. 
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FIG. 14: (color online) The interaction F for Jf = 2.60, Ji = 
— 1, J2 = —0.5. Solid lines are obtained in the second order in 
S 1-1 . The broken lines are in the second order in (Jf) -1 . The 
curves are obtained respectively for S = 00 (black), 1 (red), 1/2 
(blue) beginning at the top. The right vertical line (red) is the phase 
boundary for S = 1 where the gap closes in the second order in S~ 
at h = 0. The left vertical line (blue) is for S = 1 /2. The shift of the 
boundary near Jf = —0.5 which determines the stability of the SS 
is extremely small and is within the error of 0{S-' 2 ) orO((Jf)" 2 ). 



Finally, we comment on the shift of the boundary which de- 
termines the stability of the SS given by the MFT. As is seen 
in FigJT4j the quantum effect to T is very little on this bound- 
ary and the shift is within the error of approximation for both 
approximation (0(S~ 2 ) or 0(( Jf ) -2 )). We found that this is 
also the case for the SS2 which mainly appears for J% > 
unless the parameter set sits near the phase boundary of solid- 
SF at half filling. Therefore, if the energy of the solid phase 
is sufficiently less than the superfluid phase at half filling, the 
boundary determining the stability of the SS phase given by 
the MFT may not be affected by the quantum fluctuation. 




FIG. 13: (color online) The interaction F for Jf = 1.59, Ji = 
— 1, J2 = —0.5. Solid lines are obtained in the second order in 
The broken lines are in the second order in ( Jf The curves 
are obtained respectively for S = 00 (black), 1 (red), 1/2 (blue) 
beginning at the top. The left vertical line (red) is the phase boundary 
for S = 1 where the gap closes in the second order in S 1-1 at h = 0. 
The right vertical line (blue) is for S = 1/2. Near the boundary 
and for large Jf , the difference becomes large. For the large Ising- 
like anisotropy Jf , the evaluation of F in the second order in S~ 
becomes pathologic as discussed in the last part of Sec. lIVBl 



VII. SUMMARY 

By using the spin-wave- (1/5) and the Ising expansion to- 
gether with the dilute-Bose-gas technique, we studied the SS 
phase around the half-filled solid (NAF) phase. First, we in- 
troduced two kinds of magnon excitations for the two sublat- 
tices in the NAF phase. At a certain value of chemical po- 
tential (or, the external magnetic field), the gap of one of the 
Bogoliubov-transformed magnons closes; this magnon BEC 
keeps the two-sublattice NAF structure intact implying the SS 
phase. The spin configuration of the SS phase was determined 
by the minimum of the energy spectrum over the solid ground 
state. The Bogoliubov-transformed magnons can be viewed as 
vacancies or interstitials introduced in solids. If the effective 
interaction T among the condensed magnons are repulsive, we 
may expect, on physical grounds, a stable SS phase to appear. 
Therefore, the necessary condition for a second-order solid- 
SS transition is given by T > 0; if this condition is met, the 
SS phase realizes for low condensate density. 

To evaluate the excitation spectrum in the solid phase and 
the effective interaction V in a quantum-mechanical manner, 
we developed the perturbation theory in S" 1 and (Jf ) -1 in 
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TABLE II: Types of phase transitions suggested by the several meth- 
ods in this paper. 'SW and Tsing' represent the spinwave- and the 
Ising expansion discussed in Sec. |IV| and Sec.[V] respectively. The 
interactions Tsissi (i = 1, 2) and Tuss; are given in eqs.d33t. \34\ 
and l !51b . respectively. Fs2SSi and Fi2sS; are shown in Figs ll3ll4l 
If r diverges, one should not take the value literally, since even in 
that case a second-order solid- 'bound-magnon SS (BMSS)' transi- 
tion may be expected. The detailed discussion on a BMSS is given 
in Sec.lVCl 



method 


solid-SF 


solid-SS 


solid-BMSS 


SW 1st 


1st (Fsissi < 0) 


2nd (r S issi > 0) 




SW2nd 


1st (r S 2SSi < 0) 


2nd (r S2S Si > 0) 




Ising 1st 


1st (rnss, < o) 


2nd (Tnss, > 0) 


2nd 


Ising 2nd 


1st (ri2s Si < o) 


2nd (Tiass, > 0) 


2nd 


MFT 


1st 


2nd 





Sec. lIVI andlVl respectively. 

The first-order calculation in S^ 1 yielded the same results 
as in the MFT; three types of SS phases are found around the 
half-filled solid. At the second order in S^ 1 , we found a pos- 
sibility that quantum fluctuations destabilize the NAF solid, 
which is expected to be stable from the MFT. Specifically, in 
the evaluation of T, the second-order correction in S^ 1 be- 
comes ill-behaved when the Ising-like anisotropy Jf is large 
or when the energy of the solid is almost the same as that of 
the SF at half filling. In such cases, the MFT (or, equivalently, 
the first-order perturbation in S*" 1 ) may not be reliable. 

In order to overcome this difficulty, we carried out another 
perturbation theory from the Ising limit (i.e. expansion in 
1/Jf). At the first order in 1/Jf, we used the ladder ap- 
proximation and saw that quantum fluctuations did suppress 
T, while we obtained the same result as the MFT one as far as 
the stability of the SS phase is concerned. 

For negative T (i.e. attraction), we found a possibility of a 
novel phase characterized by the bound-magnon condensate. 
However, this phase may be replaced by the SF for the param- 
eters considered in the text. In carrying out the second order 
calculation in 1/Jf, we used the ladder approximation with 
only diagrams up to 2-loop (i.e. up to the third order in 
kept in the kernel. The effect of quantum fluctuations depends 
crucially on the energies of the solid- and the SF phase, as in 
the large-5 expansion. 

When the energy of the solid is sufficiently smaller than that 
of the SF phase at half filling, the second order term had little 
effect on T in the vicinity of the MFT-boundary (see FigfPfli. 
In other words, under the above condition, we may conclude 
that quantum fluctuations only have minor effects on the sta- 
bility of the SS phase. 

On the other hand, when the energy of the solid phase is 
comparable to that of the SF phase, there exists a possibility 
that quantum fluctuations completely wash out the SS phase 
obtained in the MFT. Actually, in the vicinity of the LF point, 
where frustration due to the competition among NAF, CAF 
and SF is strong, the second-order Ising-like expansion also 



concluded divergingly large negative values of T (see Fig. [13] 
and Fig. [Till. 

In Sec. |VTJ we studied the effect of quantum fluctuations 
on the ground state at the LF point. At the second order in 
S _1 , the ground state may be given not by the SS but by the 
SF even at h = 0. The failure of the Liu-Fisher values to de- 
scribe 4 He suggests that the optimal parameters, which should 
be obtained by fully quantum treatment, may differ from the 
Liu-Fisher ones. We expect that the deviation from the LF 
point is small and that it can be handled in a perturbative fash- 
ion. On the basis of this expectation, we studied the stability 
of the SS in the vicinity of the LF point. Since the energies of 
SS and SF are comparable in this region, the MFT may not be 
reliable. Even if the shift of the fitting parameters from the LF 
point is small and the MFT guarantees the stability of the SS, 
there remains a possibility that quantum fluctuations destabi- 
lize the SS. To investigate this possibility more closely, we 
shall need such a sophisticated treatment that the renormal- 
ization of the effective interaction V due to higher order terms 
is appropriately taken into account. 

Note added - After the completion of our work, we became 
aware of a series of papers by Stoffel and Gulacsi who stud- 
ied the same model^ as ours by the Green's function theory 
with the random-phase approximation. They reached a differ- 
ent conclusion that the critical external field at the solid-SS 
transition is little affected by quantum fluctuations at the LF 
point. We suspect that the discrepancy might be attributed to 
the difference in the approximation schemes; we believe that 
our approximation is well controlled by the two small param- 
eters. 
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APPENDIX A: SUPERSOLID PHASE EMERGING FROM 
TWO TYPES OF BOSONS 

In Sec. IIV Al we have discussed the two SS phases, which 
are described by a single Boson condensate, and classified 
them by A ( |29l . As has been mentioned there, however, in 
the case of A = 0, the spin wave dispersion takes its minima 
at both Qi = (0, 0, 0) and Q2 = (tt, it, it), and there exists 
a possibility that both kinds of bosons condense simultane- 
ously. In this appendix, we discuss this possibility within the 
first-order perturbation in S~ 1 . We shall see that a new type of 
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SS phase (SS3) appears for a certain parameter region; it has 
a 4-sublattice structure and may continue to the quarter-filled 
solid. 

As in the case of magnon BEC just below the saturation 
fiel d 22 ' 24 ' 25 , the ground-state energy density may be expanded 
in powers of the boson densities; 

E s ff 1 2 1 2 

— « const + -r Ql p Ql + -r Q2 p Q2 + T 2 p Ql PQ 2 

+ r 3PQiPQ 2 cos2(^ Ql - lfQ 2 ) - Sp a (p Ql + pq 2 ) , 

(Al) 

where 

T 2 = (V a (0; Qi, Q 2 ) + V a (Q 2 - Q i; Qi, Q 2 ) 

+ 14(0; Q 2 , Qi) + V a (Qi - Q 2 ; Q 2 , Qi))/2 , (A2) 
r 3 = V a {Q 2 ; Qi, Qi) (= V a (Q 2 ; Q 2 , Q 2 )) , 



and (a q ) = y/Np^e^ , T Q] = T sissj for j = 1, 2 (eqs.® 
and (l34li). Since the Hamiltonian ^ is not hermitian, it is not 
always true that Vq(Q 2 ; Qi, Qi) = V Q (Q 2 ; Q 2 , Q 2 ). How- 
ever, these coincides with each other when A = 0. The rela- 
tive angle ((fiQ 1 — (fiQ 2 ) takes (it/ 2) when T 3 < (> 0). 

If T Qi < or v /rQ i r Q2 < -(r 2 - |r 3 |), a magnetization 
jump occurs. Otherwise, when Min^c^ , Tq 9 ] < T 2 — |r 3 |, 
only one of the two species, which has smaller r^. condenses 
and forms the spin structure ( f20l >. If Min^Qj , TqJ > T 2 — 
1 17*3 1 , ( IA5I ) takes the minimum when 



PQi 

PQ2 



Q 2 



(T 2 



r Ql r Q2 -(r 2 -|r 3 |) 
r Ql -(r 2 -|r 3 |) 



;Sp 



(A3) 



r Ql r Q2 -(r 2 -|r 3 | 

Then, the spin configuration is given by, 

(S?) = V25(VpQTcoBhflgJ cos^ Ql 



2 <S>0 



/pcj7cos(Q 2 ' R i + ( ^Q 2 ))( 1 
(i). 



f(AS^) 



S 



(Sf) = ±V2S(y^7cosh^ sin^ Ql 



f(AS^) 

/po7sin(Q 2 • R; + <Pq 2 ))(1 + g J ),for leA 



(S?) = (S - AS™) - ( PQl cosh 2 0« 



S 



-cosh Qi ^/Pq 1 Pq 2 cos(Q 2 • Rz + <^q 2 - vqJ) , 



(A4a) 



(S*) = -V^Tsinh^ cos^ Ql (l 

(S»> = TV^Tsmh^ sin^ Ql (l 

(i 
Q 



/(A5 (1) ) , 
5 

S ' 



(S*J = -(S - AS^ 1 
for m £ B 



PQj sinh 2 6>r 



(A4b) 



where we use sinh Oq^ = and AS and f(AS) is the same as 
in eq.(|20l. By some numerical calculations, we found that this 
non-trivial SS phase with ((fQ 1 — <^q 2 ) = (SS3) is stabilized 
for a broad region of the parameter-space, mainly for J| > 0. 
For example, if Jf /| Jfl = 3 (and A = 0), the SS3 exists for 
0.2 < J$/\Jt\ < 2-0 . 

For A w 0, r 2 and T 3 may have a influence on the magne- 
tization process around the half-filled solid. For example, if 
A > 0, the system energy is given by 



-Eeff 



1. 



; r QiPQi + 2 F Q2^Q2 + r 2PQiPQ2 
r 3 PQ 1 PQ 2 cos2(«^ Ql - 9?Q 2 ) 
SpoPcn + (-Spo + A 2 )p Q2 , 



(A5) 



where A 2 = e c i(Q 2 ) — e c l(Qi) ~ O(A) > and Ts obtained 
at A = may be used approximately. If the used parame- 
ters satisfy the condition of the stability of the SS3 discussed 
above, a phase transition from SSI to SS3 occurs at: 



Spoc 



r Ql A 2 



r Ql - (r 2 - |r 3 |) 



(A6) 



Then, the densities of the condensed bosons are given by 



(rq 2 -(r 2 -|r 3 |))S/io + (r 2 -|r 3 |)A 2 rA7 , 

PQi = ^ — ^ it-i i\9 ' (A /a) 



PQ2 



r Ql r Q2 -(r 2 -|r 3 |)2 

(r Ql -(r 2 -|r 3 |))5 At0 -r Ql A 2 
r Ql r Q2 -(r 2 -|r 3 |)2 



(A7b) 



At po = Moci> ( IA7t and (l22l give the same density pq ± 2 , 
and thus a second order phase transition is implied. If A < 0, 
similarly, a second order phase transition from SS2 to SS3 
occurs at: 



Sp0c2 



Tq 2 Ai 



r Q2 - (r 2 - |r 3 |) 



(A8) 



where Ai = — A 2 . 



APPENDIX B: SOME EQUATIONS OMITTED IN THE TEXT 

1. Section lTVBl 

The additional quadratic terms in eq.d39l) emerging from 
normal order of the Bogoliubov-transformed bosons are given 
by: 
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T\(k) =e (k)(l]Tsinh 2 ^) - J 1 x (l^C 1 (q) S inli2^) + (-| J|C 2 (k) + 2J 2 L )(i £ C 2 (q) sinh 2 0,) , 
q q q 

i j« i (B1) 

T 2 (k) =^(k)(-^ S inh 2 g ) - -^C 1 (fc)(-^C 1 (q)sinh2^) , 
q q 

2. SectionlVBl 

The interaction part of Hamiltonian which contributes to the kernel of the order of ( Jf ) _1 is given by: 
#(i2)int = J7 X] {(^^(q) - J 2 L C 2 (k 2 ) - 2J 1 z sinh6» k2 _ q sinh0 k2 Ci(q) + J^C^ka) sinh6' k2 )aj ci+q aj i2 _ q a kl a k2 

q,ki,k 2 

- 2J^Ci(q)a kl+q ^ k2 _ q a kl ^k 2 + 2JfCi(q) sinh0 k2+q a ki+q a k2+q a kl /3_ k2 
+ (2J 1 z C 1 (q)sinh0 k2 _ q - J 1 X Ci(k 2 ))4 i+q 4 2 _ q a kl ^l k2 } 

(B2) 
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In Sec. IJnl the hermiticity of the Hamiltonian is assumed. In 
general, the non-hermitian DM Hamiltonian Hdm is given by 
U^ 1 HupU, where Hhp is the Holstein-Primakoff transformed 
Hamiltonian and U is the non-unitary operator, which recovers 
the hermiticity of .Hdm-~ In NAF phase, we may not consider the 
effect of U when we calculates the observables since the ground 
state (vacuum) is the eigenstate of U. If the boson condenses, U 
may shift the vacuum and we must manipulate the operator U ex- 
plicitly. However, since we consider only a dilute-gas limit and 
need the observables obtained in the NAF phase (at /j, a = CP), 



the hermiticity does not matter in the concrete discussion. 
When there exist degenerate minima at several Qs, there is a pos- 
sibility that the bosons at different Qs simultaneously condense 
and the magnetic structure may be different from the one charac- 
terized by l |20| >. See appendix A for more detail. 
The analysis near the saturation field is free from this problem; 
the fully polarized state is an exact eigenstate and high enough 
magnetic field guarantees the validity of the reference state. 



